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Abstract 

We consider the most general Dunkl shift operator L with the following properties: (i) L is of first order 
in the shift operator and involves reflections; (ii) L preserves the space of polynomials of a given degree; 
(iii) L is potentially self-adjoint. We show that under these conditions, the operator L has eigenfunctions 
which coincide with the Bannai-Ito polynomials. We construct a polynomial basis which is lower-triangular 
and two-diagonal with respect to the action of the operator L. This allows to express the BI polynomials 
explicitly. We also present an anti-commutator AW(3) algebra corresponding to this operator. From the 
representations of this algebra, we derive the structure and recurrence relations of the BI polynomials. We 
introduce new orthogonal polynomials - referred to as the complementary BI polynomials - as an alternative 
g — >■ — 1 limit of the Askey- Wilson polynomials. These complementary BI polynomials lead to a new explicit 
expression for the BI polynomials in terms of the ordinary Wilson polynomials. 
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1. Introduction 



Two new families of "classical" polynomials [32], [33] were introduced recently through limits when q goes to -1 
of the little and big q-Jacobi polynomials. The term "classical" is taken to mean that the polynomials 
are eigenfunctions of some differential or difference operator L: 



In these instances, the operator L is of first order in the derivative operator dx and contains moreover the 
reflection operator R defined by Rf{x) = f{—x); it can be identified as a first order operator of Dunkl type [10] 
written as 



with some real rational functions F{x), Fi{x), G{x). 

The eigenvalue equation (1.1) (with L as in (1.2)) has further been investigated generally in [34]. As a rule 
the polynomial eigenfunctions Pn{x) are not orthogonal. However, if F\{x) = the big and little -1 Jacobi 
polynomials are seen to be the only classes of orthogonal polynomials consisting in eigenfunctions of such first 
order Dunkl differential operators. 

The main purpose of the present paper is to study a difference analogue of the operator L and the corre- 
sponding eigenpolynomials Pn{x). We introduce a natural generalization of the operator L which (apart from 
the multiplication operator) only contains 2 nontrivial operators: the shift T+ and the reflection R. We shall 
call such operators Dunkl shift operators on the uniform grid. We flnd the conditions for such operators to 
transform any polynomial into a polynomial of the same degree. Under these conditions the eigenvalue problem 
can be posited and we can look for the corresponding eigenpolynomial solutions Pn{x). 

The main result of our study is that the polynomial eigenfunctions of such a generic Dunkl shift operator, 
coincide with the Bannai-Ito (BI) polynomials that depend on 4 parameters. 

Recall that the Bannai-Ito polynomials were flrst proposed in [3] as a g —1 limit of the q-Racah poly- 
nomials. Bannai and Ito showed that these polynomials together with the q-Racah polynomials (and their 
specializations and limiting cases) are the most general orthogonal polynomial systems satisfying the Leonard 
duality property. Bannai and Ito derived the three-term recurrence relation for their polynomials and also 
presented explicit expressions in terms of linear combinations of two hypergeometric functions 4-F3(l) [3]. 

Further developments of the theory of BI polynomials and their applications can be found in papers by 
Terwilliger [29], [30], Curtin [8], [9] and Vidunas [31]. 

Our approach is completely different: we start with a generic first order Dunkl shift operator that preserves 
the space of polynomials and then show that its polynomial eigensolutions coincide with the BI polynomials. 
This approach gives a new type of eigenvalue equation for the BI polynomials - one that involves a combination of 
the ordinary shift operator with refiections. Moreover, we give an algebraic interpretation of the BI polynomials 
in terms of a g' = — 1 version of the AW(3) algebra. This interpretation allows one to derive the structure 



LPn{x) = XnPn{x). 



(1.1) 



L = F{x){I -R) + Fi{x)dx + G{x)dxR 



(1.2) 
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and recurrence relations for the BI polynomials. We introduce also a new class of orthogonal polynomials - 
the complementary Bannai-Ito polynomials. They can be obtained from the Askey- Wilson polynomials by a 
limiting process q ^ —1 which is slightly different from the one used by Bannai and Ito in [3] to obtain their 
polynomials. In contrast to the BI polynomials, the complementary BI polynomials do not possess the Leonard 
duality property. Nevertheless, the complementary BI polynomials have a very simple expression in terms of 
the ordinary Wilson polynomials. Moreover, the complementary BI polynomials are Christoffel transforms of 
the BI polynomials. This allows one to present a new explicit expression of the BI polynomials in the form of 
a linear combination of two Wilson polynomials. 
The present paper is organized as follows. 

In the next section we introduce a class of (potentially symmetric) operators L which has the following 
properties: 

(i) L is of first order with respect to T+ and contains also the reflection operator R; 

(ii) L preserves the space of polynomials (i.e. it transforms any polynomial into a polynomial of the same 
degree). 

These properties imply the existence of eigenpolynomials Pn{x) satisfying (1.1). 

In Section 3, we determine these eigenpolynomials Pn{x) explicitly. Most useful to that end is the observation 
that there is a polynomial basis (<?!>„ (x)) in which the operator L is lower triangular with only 2 diagonals, i.e. 

L(t>n{x) = Xn(pn{x) + i/„0„_i(x). 

This property leads to an explicit formula for the polynomials P„(a;) in terms of a linear combination of two 
hypergeometric functions 4^3(1). 

In Section 4, we show that the operator L together with the operator multiplication by x, forms a g = — 1 
(i.e. an anticommutator) analogue of the ^W(3) algebra which we call the Bannai-Ito algebra. Using the 
representations of this algebra, we determine the structure relations of the BI polynomials Pn{x) and show that 
these polynomials satisfy expectedly a 3-term recurrence relation 

Pn+l{x) + hnPn{x) + U„P„_i(a;) = xPn{x), 

where the recurrence coefficients u„, 6„ are derived explicitly. This allows to identify the eigenpolynomials Pn{x) 
with the Bannai-Ito polynomials first introduced in [3]. 

In Section 5, we construct the complementary Bannai-Ito polynomials and obtain a new explicit expression 
of the BI polynomials in terms of ordinary Wilson polynomials. 

In Section 6, we demonstrate that the operator L is symmetrizable. This means that there exists a function 
ip{x) such that the operator ip{x)L is symmetric 

{ip{x)L)* = ip{x)L, 

where L* means the formal conjugate of L (constructed according to natural rules). 
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In Section 7, we introduce the Bannai-Ito grid, which is a discrete invariant set of the BI operator. The BI 
operator becomes a 3-diagonal matrix when restricted to this set. Using this observation we derive the weight 
function for the case when the BI polynomials are orthogonal on a finite set of points of the real line. 

In Section 8, it is shown how the BI and complementary BI polynomials can be obtained by limiting processes 
from the Askey- Wilson polynomials. 

In Section 9, it is indicated how the Dunkl shift operator can be derived from the Askey- Wilson difference 
operator when q ^ —1. 

In Section 10, we consider a special symmetric case of the BI polynomials. This leads to the construction of a 
new family of orthogonal polynomials (explicitly expressed in terms of dual continuous Hahn polynomials) with 
a purely continuous weight function on the whole real axis. This is the first nontrivial example of an infinite 
positive family of BI polynomials. (Recall that only finite systems of BI polynomials have been considered in 
details so far). 

In Section 11 we show that a limiting process x ^ x/h,h ^ reduces the Dunkl shift operator to a Dunkl 
differential operator L having the form (1.2) with Fi{x) = 0. This means that correspondingly, the Bannai-Ito 
polynomials tend in this limit to the big -1 Jacobi polynomials introduced in [33]. 

2. Dunkl shift operators on the uniform grid and their polynomial 
eigensolutions 

In complete analogy with the differential case [34], let us consider the most general linear operator L of first 
order with respect to the shift operator T+ which also contains the reflection operator R: 

L = Fo{x) + Fi{x)R + Gq{x)T+ + Gi{x)T+R, (2.1) 

where Fq{x), Fi{x),Go{x),Gi{x) are arbitrary functions. The shift operator is defined as usual by: 

T+ = exp(5^) 

so that T+/(a;) = f{x + 1) for any function f{x). The operator T+i? acts on functions f{x) according to 
T+Rf{x) = fi-x-l). 

We are seeking orthogonal polynomial eigensolutions of the operator L, i.e. for every n we assume that there 
exists a monic polynomial Pn{x) = x"' + 0{x"'~^) which is an eigenfunction of the operator L with eigenvalue 
A„ (1.1). In what follows we will suppose that 

A„ ^ for n = 1, 2, . . . , A„ 7^ for n 7^ m. (2.2) 

We first establish the necessary conditions for the existence of such eigensolutions. 

We will assume the operator L to be potentially self-adjoint. This leads to the condition 

^0(2;) = (2.3) 
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(see also Sect. 5 for further details). On the one hand, the formal adjoint of the operator T+ is the backward-shift 
operator (T+)* = T~ , where 

T-f{x) = f{x-l). 

On the other hand, for the adjoint of the operator T+i? we have 

[T+R]* = R*T- = RT- = T+R. 

Thus, the operator T+i? is (formally) self-adjoint, while the operator is not. Hence the symmetric (or 
potentially self-adjoint) operators of type (2.1) cannot contain a term with the operator only. 

Moreover, without loss of generality we can assume that Fq + Fi + Gi = 0. Indeed, L{1} = constant. We 
can assume that L{1} = (otherwise we can add a constant to Fo{x) to ensure this). This leads to the desired 
condition. 

We thus can present the operator L in the form 

L = F{x){I -R) + G{x){T+R - I) (2.4) 

with only two unknown functions F{x), G{x). Here / is identity operator. 

Equivalently, we can write the action of the operator L on functions f{x) as: 

Lf{x) = F{x){f{x) - f{-x)) + G{x){f{-x - 1) - f{x)). (2.5) 

Considering the action of the operator L on x and x^ we arrive at the necessary conditions 

F(.) = ^iM±^, G(.) = |l^, (2.6) 

where by gj(a;), i = 1, 2 we mean arbitrary polynomials of degree i, i.e. 

<l\{x) = ^ix + £,0, q2{x) = r]2x'^ + r]ix + r]o (2.7) 

with arbitrary coefficients • • • : ''72. 

It is easily seen that these conditions are also sufficient; namely, that the operator L defined by (2.4) with 
functions F{x),G{x) given by (2.6), (2.7), preserves the linear space of polynomials. More exactly, for any 
polynomial Q{x) of degree n we have LQ{x) = Q{x), where Q{x) is another polynomial of degree n. For 
monomials x" we find 

La;" = A„a;" + 0(a:"-^), (2.8) 

where 

if n is even 

^.-miipll if n is odd. ^^-^^ 
From this expression it is clear that necessarily 

m^O, ^ij^mN, N = 0,1,2,... (2.10) 
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(otherwise the operator L becomes degenerate: it does not preserve the degree of the polynomial). 
Assuming that conditions (2.10) are valid, we can construct the eigensolutions 

LP„(x) = A„P„(x), (2.11) 

for every n = 0, 1, 2, . . ., where Pn{x) = x" + 0(a;"~^) is a monic polynomial of degree n. 

Since r}2 7^ 0, we can always assume that 772 = 1- We can thus present the functions F{x) and G{x) as 

a^,^^i--r. + m{x-r2 + l/2)^ ^(,)^(^-pi)(^-P2) (2.12) 

with 4 arbitrary real parameters ri,r2, pi, P2- We will use this parametrization in what follows. 
The eigenvalue A„ has the expression 

{2- if n is even 
ri + r2 — pi - P2 - 2 1^ " 
From the properties of the operator L described above, it follows that the monic polynomials Pn{x) exist for 
all n = 0, 1, 2, . . . and arc the unique polynomial solutions of equation (2.11) for the given eigenvalues A„. 

The polynomials P„(x) = Pn(x:ri,r2,pi,p2) depend on 4 parameters ri,r2, pi, p2- There is however an 
obvious invariance under the action of the Klein group Z2 x Z2. Indeed, the permutations ri •<-> r2 and pi O p2 
leave the operator L unchanged and hence the polynomials P„(a;; ri,r2, Pi, P2) are invariant as well under these 
permutations. (With the identity this Z2 x Z2 group contains 4 elements.) 



3. A 2- diagonal basis and an explicit expression in terms of hyper- 
geometric functions 

Choose the following polynomial basis. For even degree, take 

<p2n{x) = {x-n + l/2)„(-a; - n + l/2)„ (3.1) 

and for odd degree, take 

hn+i{x) = {x-ri + l/2)„+i(-a; - n + l/2)„ = (x - n + l/2)(a; - n + 3/2)„(-x - n + l/2)„. (3.2) 

Here (a;)„ = x{x + l){x + 2) . . . (a; + n — 1) stands as usual for the shifted factorial (Pochhammer symbol). 
Obviously, (pn {x) is a polynomial of degree n: 

</,„(a;) = {-l)n{n-l)/2^n ^ (3.3) 

It is directly verified that the operator L is two-diagonal in this basis 

L(f)n{x) = K4>n{x) + Un(pn-l{x), (3.4) 
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where 

C §(ri+r2-f) if n is even 



{Pi ^ ^1 + 'n/2){p2 — ri+ n/2) if n is odd 
Using this striking observation, we can exphcitly construct the polynomial eigensolutions of the operator L in 
the same manner as in [33]. 

Let us expand the polynomials Pn{x) over the basis 4>n{x): 



s=0 



From the eigenvalue problem (2.11), we have the recurrence relation for the coefficients Ans'- 



Hence the coefficients Ans can easily be found in terms of A„, 



0- 



A A (An - Ao)(An - Ai) . ■ . (A„ - \s-l) , . 

1/11/2 ■■■Vs 



or in terms of the coefficient Ann'- 



= Ann ^nVn-l ■ ■ ■ ^s+l ^^ ^^ 

(An — An-l)(An — An-2) • • • (An — Ag) 

We thus have the following explicit formula for the polynomials Pn{x) 

□ / X . sr-^ (An — Ao)(A„ — Ai) . . . (A„ — As_i) ... , . 

Pn{x) =AnO > ^ ^ ^ (l>s[X). (3.9) 

U1U2 . . .Vg 

s— 

Expression (3.9) resembles Gauss' hypergeometric function and can be considered as a nontrivial generalization 
of it. 

The expansion coefficients Ans have different expressions depending on the parity of the numbers n and s. 
When n is even n = 2, 4, 6, ... we have 

A lA (-n/2),(n/2 + pi +P2 -ri -r2)s 

An,2s Ano — -TT- \ M ir, , wi /o i \~' « = 0,1,... n/2 (3.10) 

s!(l - r2 - ri)s(l/2 + Pi - ri)s(l/2 + p2 - ^1)^ 

and 

A lA c (1 - n/2), (n/2 + pi + - ri - ra). ._ni /9 i in 

An,2s+1 AnO = t,n -Tf. \ to In , \ (o let , « = U, 1, . . . n/2 - 1, (>>-ll) 

s!(l - r2 - ri)s(3/2 + pi - ri)s(3/2 + p2 - r-i)s 

where 

n. 



2(1/2 + pi-ri)(l/2 + p2-ri)" 
When n is odd n = 1, 3, 5, . . . then 

An,2s Ano = -rj- \ n /o , \ ft In , ^~' s = 0, 1, . . . (n - l)/2 (o-12) 

s!(l - r2 - ri)s(l/2 + Pi - ri)5(l/2 + p2 - ri)s 



and 



A /A - n ((1 - n)/2)si{n + 3)/2 + pi + pa - n - ra). , , . 

An,2s+1 AnO — -^n -rr-. \ to In , \ to In i T' S — I, . . . {U - I) I, {6.16) 

s!(l - r2 - ri)s(3/2 + pi - ri)s(3/2 + p2 - ri). 
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where 

(n + l)/2 + pi+ P2 - ri - r2 



Vn 



(1/2 + pi -ri)(l/2 + p2-ri)' 
It is now convenient to separate the even and odd terms in the expression for the polynomials: 

Pn{x) = ^An,2s(t>2s{x) + ^ An^2s+lhs+l{x). (3.14) 

Taking into account the explicit formulas (3.10)-(3.13) for the expansion coefficients, we finally obtain the 
following expressions in terms of hypergeometric functions: 

(i) if n is even 

Ano ^ H l-ri-r2,l/2 + pi-ri,l/2 + p2-ri ' ) 

(3.15) 

/l-n/2,n/2+l + /Ji+p2-ri-r2,x-ri+3/2,-x + ri + l/2 
V 1 -ri -r2,3/2 + /9i -ri,3/2 + p2 -ri 

(ii) if n is odd 

Pn{x) ^ /-(n-l)/2,(n+l)/2 + pi+p2-ri-r2,a;-ri + l/2,-a; + ri + l/2 \ _ 
Ano ^ l-ri-r2,l/2 + pi-ri,l/2 + p2-ri ' ) 

(3.16) 

( ^1/0^ p /^-(ri-l)/2,(n + 3)/2 + pi+p2-ri-r2,x-ri+3/2,-x + ri + l/2 \ 

rin {x-n + 1/2) 4i^3 , „ ,„ , „ ,^ , ; 1 • 

\ 1 - ri - r2,3/2 + Pi - ri,3/2 + p2 - ''I / 

Note that we have terminating hypergeometric functions (i.e. only a finite number of terms in the hyperge- 
ometric series appear). Hence, as is easily seen, the expressions in the rhs of (3.15) and (3.16) are polynomials 
of degree n in the argument x. The coefficient A„o (ensuring that the polynomials P„(a;) are monic) can be 
found from the coefficient in front of the leading term x". 

Remark . At first sight one might be tempted to identify the hypergeometric functions 4-^3(1) with the 
Wilson (Racah) polynomials which are known to be expressible in terms of truncated 4^3(1) functions [17]. 
This is not so, however, because the balance is not right. 
Recall that the generalized hypergeometric function 

ai, 02, . . . Op _ 
61,62, ' 

is said to be fc-balanced if 01+02 + . . .+ap = A;+6i+62 + . . .+69. The Racah polynomials are expressed in terms of 
-1 balanced hypergeometric functions 4.^3(1) [17]. In (3.15), (3.16), as is easily seen, all hypergeometric functions 
4^3(1) are instead zero-balanced. In what follows we will give another explicit formula for the polynomials Pn(x) 
which will relate them to "true" Wilson polynomials. 

Remark that the basis ^„(a;) is important in Terwilliger's approach to Leonard duality [29], [30]. Similar 
bases exist for -1 Jacobi polynomials as well [33]. 

Equivalently, the basis 4>n{x) can be presented as 



^„(x) = (-l)"("-i)/2[](x-ai), 
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where 



„„ = (-ir|n + M)-i. 



This means that the (j)n{x) form a special case of the Newtonian interpolation polynomial basis = 
nr=i('^ — c^i), where are fixed interpolation nodes [13]. In the theory of orthogonal polynomials sueh bases 
were first proposed by Geronimus [11]. All known explicit examples of orthogonal polynomials (including the 
Askey- Wilson polynomials at the top level) have simple expansion coefficients with respect to a Newtonian basis 
Cln{x) with appropriately chosen nodes [29], [30]. We thus see that the BI polynomials satisfy this property 
as well. 

4. Bannai-Ito algebra, structure relations and recurrence coefficients 
for the BI polynomials 

In this section we derive the 3-term recurrence relations as well as the structure relations for the Bl-polynomials. 
The main tool will be the Bannai-Ito algebra (Bl-algebra for brevity) which is a special case of the AW(3) 
algebra introduced in [35]. For the Askey- Wilson polynomials the structure relations could be derived in a 
similar manner using representations of the AW(3) algebra [18] , [15] . For a study of the structure relations of 
the -1 Jacobi polynomials see [28]. 
Consider the operators 

X = 2L + K, Y = 2x + l/2, (4.1) 
where by x we mean the operator multiplication by x and where 

K = pi + p2 - n - r2 + 1/2. 

We define also the third operator 

^_ nr2 piP2 {x-pi){x-p2) (x-ri + l/2)(x-r2 + l/2) ^^ 

a; -I- 1/2 a; x x + 1/2 ' ^ ' ' 

It is easily verified that the operator Z transforms any polynomial of degree n into a polynomial of degree n + 1 
and that 

Zx" = 2(-l)"+ia;"+i + 0(a;"). (4.3) 

The operators X, Y, Z are seen to have the following Jordan or anticommutator products which are taken 
to be the defining relations of the BI algebra: 

{X,Y} = Z + iJs, {Z,Y} = X+cou {X,Z} = Y + uj2, (4.4) 

where 

wi = 4(pip2 + r-ir2), co2 = 2{pI+ pj-rl-rl), C03 = A{pip2 - rir2)- 
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(Note that the first relation (4.4) can be considered as a definition of the operator Z). 

In the simplest case uji = U2 = oj^ = 0, the BI algebra was considered in [1], [12], [25], [26] as an anticom- 
mutator analogue of the rotation algebra so(3). A contraction of (4.4) has also been used in [28] to provide an 
algebraic account of structural properties of the big -1 Jacobi polynomials. 

The Casimir operator Q commuting with all generators X, Y, Z of the BI algebra has a simple form: 

Q = + Y'^ + Z"^. (4.5) 

In the given realization of the operators X, Y, Z in terms of difference-reflection operators, this Casimir operator 
takes the value 

Q = 2{pj + pl+rl+rl)-l/4. (4.6) 

The BI algebra can be exploited to find explicitly the 3-term recurrence relation of the polynomials Pn{x). 
We first introduce two important operators J+ and J_ by the formulas 

J^ = iY + Z){X-l/2)-^^^±^ (4.7) 



and 



J- = {Y-Z){X + l/2) + '^^^^. (4.8) 



From the commutation relations (4.4), we find that the operators J± satisfy the (anti) commutation relations 

{X,J+} = J+, {X,J_} = -J_. (4.9) 

Moreover, from (4.9) it is seen that both and commute with the operator X: 

[X,Jl] = [X,Jl]=0. (4.10) 

The operator J+ has the following property: it annihilates any constant J+{1} = and transforms any monomial 
x" of even degree n into a polynomial of the same degree n (providing that ri + r2 ^ n/2): 

J+x" =n(n-2ri -2r2)a;"-|-0(a;"-^), n = 0,2,4,..., (4.11) 

while any monomial x" of odd degree is transformed into a polynomial of even degree n + 1 (providing that 

ri + r2 - pi - P2 ^ n + I) 

J+a;" = 4(ri+r2-pi-p2-n-l)a;"+^+0(a;"), n = 1, 3, 5, . . . . (4.12) 

Quite similarly, the operator J_ transforms any monomial of even degree n into a polynomial of odd 
degree n + 1: 

J_x" = 4(pi+p2-ri-r2+n + l)x"+i+0(2;"), n = 0,2,4,..., (4.13) 
while the monomial x" of odd degree is transformed to a polynomial of the same degree n: 

J_x" = (4(pip2-rir2)-n2 + 2n(ri+r2))x" + 0(a;"-i), n = 1, 3, 5, . . . . (4.14) 
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A careful analysis (where the leading and next to leading monomials are considered), shows that the operator J+ 
transforms any polynomial of even degree n into a polynomial of degree n or less, while it maps any polynomial 
of odd degree into a polynomial of exact degree n + 1; similarly, it is seen that the operator J_ transforms any 
polynomial of even degree n into a polynomial of exact degree n + 1, while it maps any polynomial of odd degree 
into a polynomial of degree n or less. 

Now, let tpn{x) be any eigenfunction of the operator X: 

where 

M„ = 2A„ + K=(-l)"(n + K). (4.15) 

It then follows from relations (4.9) that the function ijjn{x) = J+ipn will again be an eigenfunction of the 
operator X corresponding to the eigenvalue jln = ^ — l^n- It is seen that fin = fin-i if n is even and that 
An = Mn+i if ^ is odd. 

A similar property is found for the operator J_. In this case the function J-ipn will be an eigenfunction of 
the operator X with eigenvalue —1 — /i„ which is for even n and /i„_i for odd n. 

We have already established that the operators J± transform polynomials into polynomials. Hence from the 
above observations, 

T ur\ j a^n^Pn-i{x), if n even, 

J+Pn{x) = <^ (4.16) 
aii'Pn+i{x), if n odd 



and similarly. 



/3i°'p„+i(a;), if n even 
l3i^^Pn-iix), if n odd 



The coefficients in the rhs of these formulas can be found by comparison of highest-order terms: 

^(0) ^ 2n{p, +p, + n/2)in +r,- n/2)(r, + r. - - p. - n/2) , = 4(,^ + _ _ _ , _ 1) 

ri+r2- pi- P2-n 

and 

= -Ain+r2-p,-p2-n-l) = -ai!\ = -n+ n/2){p, - + n/2)(p, - r, + n/2)(p, - + n/2) 

pi+ p2-ri-r2+n 

Note that formulas (4.16) and (4.17) show that the operators J± are block-diagonal in the basis of the 
polynomials Pn{x) with each block a 2 x 2 matrix. 

Prom these formulas it follows that the operators J|. have the polynomials Pn{x) as eigenfunctions 

r2p^^ /«n^4-i-Pn(a;), if n even 

J+Pn{x) = i (0) (1)„ , , (4-18) 



a^lian'Pnix), if n odd 



and similarly 

t2 



72 ^ jl3^n^Pl^liPn{x), if n even 

JiPn{x) = \ (4.19) 

I /5„_i/3n '-P„(a;), if n odd 



11 



Now, from the commutation relations (4.4) we find 

Jl = {Y + Zf{-X^ + X - 1/4) + {uj2 + W3)V4- (4.20) 
Using formulas (4.5) and (4.4) we also see that 

{Y + Zf=Y^ + Z'^ + {Y,Z} = Q-X^ + X + uji. 

We thus have 

Jl = iQ-X^+X + uJi){-X^ + X- 1/4) + (^2 + iosf/i, (4.21) 

where Q is the Casimir operator (4.5). Taking into account that in our realization of the BI algebra the Casimir 
operator takes the value (4.6), we arrive at the expression 

Jl = {{X + P1+P2- 1/2)2 _ ^ ^^^2^ - pi - P2 - 1/2)2 _ ^ ^2^2) _ (4 22) 

Quite similarly we obtain 

J! = {{X + P2-P1 + 1/2)2 _ _ ^^)2) + Pi - p2 + 1/2)2 _ _ ^^-,2) _ (4 23) 

We see that both operators J\ and J2 (which commute with the operator X) are expressible as 4-th degree 
polynomials in X. 

Starting from relations (4.16) and (4.17) one can now derive many useful relations for the polynomials Pn{x). 
As a first example we present the structure relations of the polynomials Pn{x). 
Define the operator Un \ n = 0, 1, 2, . . . as 

(i)_rj+, if n even 
^" ~ 1 J_, if n odd • ^^■'^^> 



We have in this case 



if n even 



\ Pn \ II n odd 

Similarly, define the operator Un \ n = 0, 1, 2, . . . as 

(2) fJ-, if n even 

^" 1 J+, if n odd • ^ ^ 



Then we have 



C/<2)p„(,) = e(2)p„+i(.), e« = (^V " (4.27) 

Q!„ ' , it n odd 



It is seen that the operator uli"^ plays the role of the lowering operator, while the operator Un'' serves as the 
raising operator for the polynomials P„(.t). 

We are now ready to derive, as a second example, the 3-term recurrence relations of the polynomials Pn{x). 

To do this we consider the operator 

V = J+{X + 1/2) + J_{X - 1/2) = 2y(X2 - 1/4) - wsX - W2/2. (4.28) 
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By construction, this operator is 2-diagonal in the basis of the polynomials Pn{x): 

, , f(A*„ + l/2)a<°V„_i(x) + (M„-l/2)/3i"'P„+i(x), if n even 
VPJx) = < (u tu , (4.29) 

\ i^in-l/2)^3i^>Pn-l{x) + i^in + l/2)a)^'Pn+l{x) if Tl odd ' 

where /i„ is the eigenvalue of the operator X given in (4.15). 

Moreover, from (4.1), the operator Y coincides (up to an affine transformation) with the multiplication by 
X, that is, in the polynomial basis Pn{x), we have YPn{x) = {2x + l/2)P„(a;). Hence, by the second equality in 
(4.28), 

VPnix) = HpI - l/4)(4x + 1) - Uslln " L02/2)P„{x). (4.30) 

Comparing (4.29) and (4.30), we arrive at the 3-term recurrence relation for the polynomials Pn{x): 

xPn{x) = Pn+l{x) + hnPn{x) + UnPn-l{x), (4.31) 

where the coefficients are 

^. = -1/4+^^^^^^ (4.32) 

and 

jjf^, if n even 
un = ^ . (4.33) 

{ ife' if " odd 

We can also present the expressions for the coefficients u„ in a more detailed form. For even n = 0, 2, 4, ... we 
have 

_ n (n + 2 pi + 2 Pa) (n - 2 n - 2 rz) (n + 2 pi + 2 p2 - 2 n - 2 ra) 
16 (n + pi + p2 - ri - ra) 

For odd n = 1, 3, 5, ... we have 

(n + 2pi -2ri)(n + 2pi -2r2)(n + 2p2 -2ri)(n + 2p2 -2r2) 
16 (n + pi+p2-ri -r2) 

We can now identify these recurrence coefficients with those of the Bannai-Ito polynomials. 

Recall that the monic Bannai-Ito polynomials can be defined through the 3-term recurrence relation (4.31), 

with the recurrence coefficients [3], [33] 

Un = An-lC'n, bn = 9o - A„ - C'n, (4.36) 

where 

_ rHMn±ilZ^i^+i±B), n even 

- 1 ^Mn+r-g(n+i-r3) „ „dd ^^-^^^ 

and 

Cn=\ '"7w M . (4.38) 

I _ 2Hn-r,-s'){n-r2-s') „ odd ^ ' 

In these formulas the parameters n , r2 , T3 , s* are the 4 essential parameters of the Bannai-Ito polynomials while 
the parameters h and 9o define an affine transformation x — > {x + 9o)h~^ of the argument of the polynomials 

Pn{x). 
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An easy analysis shows that the coefficients (4.32) and (4.33) coincide with the recurrence coefficients (4.36) 
under the identifications 



Oo = pi,h= 1/4, n = 2(pi - n), T2 = 2(pi - ra), T3 = -2(pi + pa), s* = 2(ri + rs - pi - pa)- (4.39) 

Note that apart from the solution (4.39), there exist 3 other solutions which can be obtained from (4.39) by the 
permutations ri O ra, pi ^ P2 ■ This corresponds to the previously mentioned symmetry of the BI polynomials 
with respect to the Z2 x Z2 group. 

We thus see that the orthogonal polynomials Pn{x) coincide with the general Bannai-Ito polynomials. 

For the reader's convenience we record the explicit formulas for the recurrence coefficients: 

Pn+l{x) + (pi -An- Cn)Pn{x) + An-lCnPn-l{x) = xP„{x), (4.40) 



where 

(ra+l+2 pi-2ri)(w+l+2pi-2r2) 



-f7 — n 1 — i — s , n even 

4:{n+l—ri—r2+pi+P2) ' 
(TO+l-2ri-2r2+2)Ji+2p2)(n+l+2pi+2p2) ^ 
4(n+l— ri-r2+pi+P2) ' 



(4.41) 



and 

2(n — 2 ri —2 r2) 



r -i 4(n-ri-r2+pi+P2)' " ^'^^'^ . - 4„x 

" S . (n-2.2+2p2)(n-2..+2p2) „ ^^d. " ^ ^ ^ 

4(n-ri-r2+pi+)92) ' 

Consider now the question of the positive definiteness of the polynomials P„(x) obeying a recurrence relation 

of the form (4.31). Recall that polynomials Pn{x) are called positive definite if u„ > for all n = 1,2, This 

property is equivalent to the existence of a positive measure dp,{x) on the real axis such that the orthogonality 
property reads 



6 

Pn{x)Pm{x)dp{x) = h„, S„m (4.43) 



where ft„ = uiU2 ■ ■ ■ Un, and where the integration limits a, b can be either finite or infinite. 

In our case it is seen that u„ < for sufficiently large n. It is hence impossible to ensure positivity for all 
n = 1, 2, Nevertheless it is possible to obtain a finite set of positive definite orthogonal polynomials. 

If Uj > for i = 1, 2, . . . , A'^ and un+i = 0, it is well known that we have a finite system of orthogonal 
polynomials Po{x),Pi {x),. . . , Pn{x) satisfying the discrete orthogonality relation 

N 

^ WsPn{Xs)Pm{Xs) = K 5„m, K = U1U2 -..Un, (4.44) 
s=0 

where Xs, s = 0,1, . . . ,N are the simple roots of the polynomial Pn+i{x). (The fact that the roots Xs are simple 
is guaranteed by the condition Un > 0, n = 1,2, ... ,N [6].) 
The discrete weights Wg are then given by the formula [6] 

P'n+i{^s) PN{Xs)Pj^_^_i{Xs) ^ ^ 
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where PA (x) are associative polynomials satisfying the recurrence relation 

Pi%ix) + bn+lPi'\x) + Un+^Pi'lM = XP^^\X) 

with initial conditions Pq^^ = 1, P'^^\x) = x — hi. 

We thus have positive definite polynomials Pn{x) which are orthogonal on the finite set of distinct points 

Xo,Xi,. . .Xn- 

Let us consider when such a situation takes place for our polynomials. 

First, assume that N is even. We have then from (4.35) that the condition un+i = is equivalent to one of 
4 possible conditions 

2{ri- pk) = N + 1, i,k = 1,2. 
We restrict ourselves with the condition 

2(r2 -P2)=N + 1. (4.46) 

Then it is sufficient to take 

N N-1 
r2=ri+e+—, pi=ri+e + d+ P2=n- 1/2 + e, (4.47) 

where ri, e, rf, are arbitrary positive parameters. Assuming (4.47), we can present the recurrence coefficients u„ 

in the following form. For even n = 2, 4, . . . , A'' 

_ n(n + 4ri - 2 + 4e + 2d + A') (-n + 4ri + 2 e + A') (n-2 + 2e + 2rf) 

16 (n - 1 + e + rf)^ 

for odd n= 1,3,..., A/'- 1 

_{n-l + 2e + 2d + N){n-l + 2d){n-l + 2e){-n+l + N) 

16 (n - 1 + e + d)^ 

From these expressions it is scon indeed that the positivity condition u„ > is satisfied for n = 1,2, . . . , N . 

Consider now the case when A > 1 is odd. Then from (4.34), we see that the condition mjv+i = is 
equivalent to one of three conditions: 

(i) pi+P2 = AN+l)/2, 

(ii) n + r2 = (A + l)/2, 

(iii) pi+ P2-ri-r2 = -{N + l)/2. 

Condition (iii) leads to singular u„ for n = {N + l)/2. Hence only conditions (i) and (ii) are admissible. 

Consider condition (ii) for dcfinitcncss. 

It is convenient to introduce the parametrization 

ri = (l-C-W2, r2 = (?7 + C + ^)/2, pi = (r? - OA P2 = {C - V - 2^ - N + l)/2, 

where ^, rj, ( are positive parameters with the restriction ^ > r/. Condition (ii) holds automatically and the 
recurrence coefficients become 

_n{-n + N -l + 2^){-n + N +l){-n + 2N + 2(,) 
16 (-n + AT + O 
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for even n and 

{n + 27]-l){-n + 2C + N){-n-2C + N + 2^){-n + 27j + 2N -1 + 2^) 



Un = 



16 (-n + 7V + 0^ 



for odd n. Again it is visible that > when n = 1, 2, . . . , A'' and un+i = 0. 

5. Complementary Bannai-Ito polynomials and an alternative ex- 
pression in terms of Racah polynomials 

Let us start with the following simple lemma 

Lemma 1 Let Pn{x) be monic orthogonal polynomials satisfying the recurrence relation 

Pn+l{x) + {9 -An- Cn)Pn{x) + C„ A„_iP„_i (x) = xPn{x), (5.1) 

with the standard initial conditions 

Po = 1, Pi(a;) =a;-^ + ^o. 

Assume that the real coefficients An, Cn are such that An-i > 0, C„ > 0, n = 1,2 . . . and that Cq = 0. Take 9 
to be an arbitrary real parameter. 
Define the new polynomials 

Pn{x) = (5.2) 
Then the monic polynomials Pn{x) are orthogonal and satisfy the recurrence relation 

Pn+l{x) + {9- An- Cn+l)Pn{x) + C„A„P„_i(x) = xPn{x). (5.3) 

The inverse transformation from the polynomials Pn{x) to the polynomials Pn{x) is given by the formula 

Pn{x) = Pn{x) - CnPn-l{x). (5.4) 

To the best of our knowledge this Lemma is due to Karlin and McGregor [16]. 

Another interpretation of this Lemma consists in the observation that formulas (5.2) and (5.4) are equivalent 
to the Christoffel and Geronimus transformations of orthogonal polynomials which in turn are special cases of 
the more general rational spectral (or Darboux) transformations [36], [4]. 

If the polynomials Pn{x) are orthogonal with respect to a linear functional a: 

{a,Pn{x)Pm{x)) =0, n + m, (5.5) 
then the polynomials Pnix) are orthogonal with respect to the functional (a; — ^)cr, i.e. 

l,a,{x-9)Pn{x)Pn,{x)) = ^, ni^m. (5.6) 
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Note also that 



which can easily be verified directly from (5.1). 

The next Lemma will be useful in the identification of polynomial systems with known families of orthogonal 
polynomials. 

Lemma 2 Let Pn{x) and Pn{x) be two systems of orthogonal polynomials defined as in the previous Lemma. 
Assume moreover that = > 0. 

Define the following monic polynomials 

S2n{x) = Pn{x^), S^n+li^) = {x - x)Pn{x'). (5.8) 

The polynomials Sn{x) are orthogonal and satisfy the recurrence relation 

Sn+l + {-lTxSn{x) + VnSn-l{x) = xSn{x), (5.9) 

where 

V2n = -Cn, V2n+1 = -^n- (5.10) 

In particular, when x = 0,. the polynomials Sn{x) become symmetric, i.e. S'„(— :r) = ( — l)"S'„(a;). 

Conversely, assume that polynomials Sn{x) satisfy the recurrence relation (5.9) with some coefficients u„ 
and a real constant x- Then the polynomials Pn{x) and Pn{x) defined by (5.8), are orthogonal and obey the 



recurrence 



relations (5.1) and (5.3) with 9 = x^ . 



This Lemma is due to Chihara [5]. For further development and applications of this Lemma see, e.g. [23], [33]. 

We shall call the polynomials Pn{x), the companion polynomials with respect to Pn{x) (sometimes the 
polynomials Pn{x) are referred to as the kernel polynomials [6]). 

The importance of Lemma 1 lies in the fact that almost all known classical orthogonal polynomials in 
the Askey table [17] admit a representation of their recurrence relation in the form (5.1) with simple explicit 
coefficients An,Cn. The corresponding companion polynomials Pn{x) satisfy similar recurrence relation (5.3) 
and belong to the same class of orthogonal polynomials (with shifted parameters). For the polynomials at the 
top of the Askey hierarchy - the Askey- Wilson (or g-Racah) polynomials, this was first observed by Chihara in 
[7]. Further development of this subject in terms of Darboux transformations can be found in [27]. 

We see that the Bannai-Ito polynomials admit as well this representation with the coefficients A„, C„ given 
by (4.41) and (4.42). 

We can hence introduce their companion polynomials, which we denote by Wn{x), as follows 

X-pl 
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with An given by (4.41). The polynomials W„(a;) satisfy the recurrence relation (5.3). Specifically, this recur- 
rence relation reads 

Wn+i{x) + (-l)>2W„(a;) + VnWn{x) = xWn{x), (5.12) 

with 

n(n + pi - ri + l/2)(n + pi - r2 + l/2)(n - n - r2) 



V2n 



{2n+l+g){2n + g) 



{n + g + l){n + pi+p2 + l){n + p2-ri + l/2){n + p2-r2 + l/2) 

V2n+i — 77, — r~i~, — — ; — r^vv ' [0.16 

{2n + 1 + g){2n + g + 2) 

where we denote 

g = pi + P2 - n -r2. 

It is seen that the companion polynomials Wn{x) satisfy a recurrence relation of the type (5.9) with x = P2- 
Hence a pair of polynomials Un{x), Vn{x) can be constructed via 

Un{x'') = W2nix), K(a;^) = n = 0,1,2,.... (5.14) 

X — p2 

(It is seen that W2n{x) are polynomials depending only on x^ and that the polynomials W2n+i{x) all have in 
common the factor x — p2)- 

From the above Lemmas, it follows that the polynomials Un{y), Vn{y) satisfy the following system of relations 

Un{y) = Vn{y) + V2nVn-l{y), (t/ - pl)Vn{y) = C/„+l(t/) + V2n+lUn{y) (5.15) 

from which the three-term recurrence relations for the polynomials Un{y) and Vn{y) can be derived: 

Un+l{y) + {pI + V2n + V2n+l)Un{y) + V2nV2n-lUn-liy) = yUn{y) (5.16) 

and 

Vn+l{y) + {pI + V2n+2 + V2n+l)V„{y) + V2„V2n+lVn-l{y) = yVn{y). (5-17) 

From the explicit expressions (5.13) of the coefficients w„, it is possible to conclude that the polynomials 
Un{x^), Vnix"^) are Wilson polynomials [17]. Omitting technical details we have 

[7„(^2) = ,F3 f -n,n + g + l,p2 + x,p2-x \ 

^ ^ " * V/0i+p2 + l,p2-ri + l/2,p2-r2 + l/2' ; ^ ^ 

and 

/ 2n (2) TP f -n,n + g + 2,p2 + l+x,p2 + l- X \ 

Vn(X ) = kI, ' 4J:<3 { „ „. „ ,„ ; 1 , (5.19) 

^ ^ " Vpi+P2 + 2,p2-ri+3/2,p2-r2 + 3/2' J' ^ ' 

where the normalization coefficients (making the polynomials monic) are 

^(1) _ {1+ Pl+ P2):,{P2 - n + l/'2)nXp2 - '1 + (2) _ (2 + + P2)„{p2 - M + 3/2),,, (p2 - M + 3/2)„ 



{n + g + l)n ' " {n + g + 2)n 

(Note that the polynomials Vn{x^) are obtained from the polynomials Unix"^) by simply shifting the param- 
eter p2 by 1: P2 P2 + 1-) 



18 



We thus have an exphcit expression for the polynomials W„(a;) (which are Christoffel transforms of the 
Bannai-Ito polynomials) in terms of the Wilson polynomials: 

W2n{x) = Un{x^), W2n+l{x) = {x - P2)Vn{x^), (5.20) 

where the polynomials Un{x'^),Vn{x'^) are given by (5.18), (5.19). 
Hence, by (5.4), we get for the Bannai-Ito polynomials 

Pn{x) = Wn{x) - CnWn-l{x), (5.21) 

where C„ is given by (4.42). Formula (5.21) can be considered as an alternative explicit expression of the 
Bannai-Ito polynomials. Like in (3.15) and (3.16), we have a presentation of the Bannai-Ito polynomials Pn{x) 
as a linear combination of two hypergeometric polynomials 4F3(1). In (5.21) however, the hypergcometric 
functions are -1 - balanced and are thus bona fide Wilson polynomials. Note that the representation (3.15), 
(3.16) is equivalent to the one obtained by Bannai and Ito for their polynomials by using a direct limit process 
from the q-Racah polynomials [3]. The representation (5.21) in terms of the Wilson polynomials seems to be 
new. When A'' — >■ oo, the Bannai-Ito polynomials tend to the big -1 Jacobi polynomials [33]. Formula (5.21) 
then becomes the formula that gives the expression of the big -1 Jacobi polynomials in terms of the ordinary 
Jacobi polynomials [33]. 

We shall refer to the polynomials Wn{x) as the complementary Bannai-Ito (CBI) polynomials. They are 
Christoffel transforms of the Bannai-Ito polynomials. They have a simpler structure (e.g. for every n the 
polynomial W„(x) is expressed in terms of only one hypergeometric function 4F3(1)). However, in contrast to 
the Bannai-Ito polynomials, the CBI polynomials do not satisfy the Leonard duality property. In other words, 
they are not eigenpolynomials of a Dunkl shift operator of the form (2.4). Indeed, the CBI polynomials have 
recurrence coefficients that do not belong to the BI class, and hence they do not arise in the classification of all 
orthogonal polynomials with the Leonard duality property realized by Bannai and Ito [3] . 

6. Symmetry factor for the Bannai-Ito operator 

Let us return to the Bannai-Ito operator (2.4), where F{x) and G{x) are given by (2.12). Consider the formal 
Lagrange adjoint L* of L: 

L* = F{x)I - R*F{x) + (T+RYGix) - G{x)I. 

We will assiime that R* = R and T+* = T' (these are natural and conventional conditions for difference 
operators on the real line). 

There are in addition, the obvious relations 

T-R = RT+, T+Rf{x) = f{-x - 1)T+R, 
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where f{x) is an arbitrary function of the real variable x. Prom these relations it follows that the operator T~^R 
is formally self- adjoint: {T+R)* = T+R. 
We then obtain 

L* = F{x)I - F{-x)R + G{-x - 1)T+R - G{x)I. 
It is visible that the Bannai-Ito operator is not symmetric, i.e. L* ^ L for any choice of the parameters 

Nevertheless, it is possible to find a symmetry factor <p{x) for the operator L such that 

{^{x)Lr = ^{x)L. (6.1) 

An operator L is called symmctrizablc if there exists a real- valued functions (p{x) ensuring that property (6.1) 
holds. The similar approach for differential operators and their polynomial solutions is familiar [22]. 

Let us show that the Bannai-Ito operator is symmetrizable for generic values of the parameters ri, r2, pi, P2- 

We have 

^p{x)L = ip{x)F{x) - ip{x)F{x)R + ip{x)G{x)T+R - ip{x)G{x) (6.2) 

and 

{<fi{x)L)* = <fi{x)F{x) - (p{-x)F{-x)R + <p{-x - l)G{-x - 1)T+R - ip{x)G{x). (6.3) 
Comparing (6.3) and (6.2), we see that condition (6.1) is valid iff the following two conditions 

ip{-x)F{-x) = ip{x)F{x), ip{-x - l)G(-x - 1) = ip{x)G{x) (6.4) 

are fulfilled. 

Clearly, these conditions are equivalent to 

ip{x)F{x) = Eo{x), ip{x)G{x) = Ei{x + 1/2), (6.5) 

where Eq{x),Ei{x) are even functions, i.e. Eq(—x) — Eq(x), Ei{—x) = Ei{x). From (6.5) we have 

F(x) E^{x) 



as well as 



and 



G{x) Ei{x + l/2) 
F{-x) Eo{x) 



G{-x) Ei{x-l/2) 
F{-x-l) Eo{x + l) 



G{-x-l) Ei{x + l/2)' 
Prom (6.6) and (6.7), we get an equation for the unknown function Ei{x) 

Ei{x + l/2) _ F{-x)G{x) 
Ei{x-l/2) ~ F{x)G{-x) ■ 



(6.6) 

(6.7) 
(6.8) 

(6.9) 
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Similarly, from (6.6) and (6.8), we find an equation for the unknown function Eo{x): 

Eo{x + l) _ F{-x-l)G{x) 



(6.10) 



Eo{x) F{x)G{-x-l)' 

The general solution of equations (6.10) and (6.9) can easily be obtained from the expressions for the functions 
F{x) and G{x): 

aoix) r(l/2 - n + x)T{l/2 - n - a;)r(l + pi + x)T{l + pi - x) 
"^""^ X r{l/2 + r2+x)r{l/2 + r2-x)r{-p2+x)r{-p2-x) ' ^ ^ 

a.jx) r(l - n + a:)r(l - n - a:)r(l/2 + Pi + x)r(l/2 + Pi - a:) 
^ a; r(r2 + x)r(r2 - x)r(l/2 - p2 + .T)r(l/2 - p2 - x) ' ^ ^ 

where (7o(a;) and cri(a;) are arbitrary functions with period 1: ao.i{x + l) = ao^i{x). The fact that both functions 
Eo{x) and Ei{x) must be even, requires that both ao{x) and ai{x) be odd: 

<^o,i{-x) = -ao^i{x). 

From condition (6.6) it follows that ai{x) = ao{x). Then, from (6.5), we find the function (p{x) 

, , T{pi- x) r jx - n + 1/2) r {~x - n + 1/2) r (x + 1 + pi) 

ip{x} = -2ao{x) — — — — — ■ , . rn^ TFT (^-^^^ 

r (.-E + 1 - P2) r [x + r2 + 1/2) r (r2 + 1/2 - xjT (-X - P2) 

Conversely, assume that the function (f{x) is given by (6.13). It can then be easily verified that the function 
(p{x) is a solution of conditions (6.4) iff the function ao{x) satisfies the two conditions: 

(i) be 1-pcriodic, i.e. aQ{x + 1) = ctq^x); 

(ii) be odd, i.e. (7q{~x) = —ao{x). 

Hence wc succeeded in constructing the desired symmetry factor 95 (.t). 

It should be noted that the function (p{x) can be presented in several equivalent forms if one exploits the 
classical formula 

T{z)r{i-z) = ^. 

simrz 

For example, one can present ip{x) in the form 

./ X _ /X r(pi ~ x)r{pi + 1 + a:)r(p2 - a;)r(l + p2 + x) 

nx) ''^^'') r{r2 + l/2 + x)T{r2 + l/2-x)T{n + l/2 + x)r{n + l/2-xy ^ ' 

where the function a2{x) has the same properties (i) and (ii) as ao{x). 

Let us introduce the operator M = ip{x)L, where <p(x) is given by (6.13). This operator will thus be formally 
self-adjoint, or symmetric: M* = M. 



7. Bannai-Ito grid and weight function for the Bannai-Ito polyno- 
mials 

One of the most important property of Dunkl shift operators like L or M is that there is a discrete set of points 
on the real line which is invariant under the action of these operators. 
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Indeed, let us introduce the Bannai-Ito (BI) grid of the first type Xg as a discrete set of real points defined 
by the two conditions 

-X2s = X2s-l, -I - X2s = X2s+\, S = 0, ±1, ±2, . . . . (7.1) 

It is easy to verify that the general solution for the grid Xg satisfying conditions (7.1), is 

{a + I if s is even, 
-a - ^ if s IS odd, 

where a is an arbitrary real parameter. 
Equivalently, 

a;, = -1/4+ (-1)^(1/4 + o + s/2). (7.3) 
Using properties (7.1), we can conclude that for any Dunkl shift operator 

H = A{x)R + B{x)T+R + C{x) (7.4) 

with arbitrary functions A{x),B{x),C{x), we have for any function f{x) 

Hf{Xs) = ^sf (Xs+l) + VsfiXs) + Csf{Xs-l), (7.5) 

where 

f = /^(^«) if sis even, 

\ A{xs) if sis odd, ^ ■ ' 

^ (A{xs) if s is even, 

\ B{xs) if s is odd, ^ ' 

% = C{xs). (7.8) 
In other words, the Dunkl shift operator becomes a 3-diagonal matrix (either finite or infinite) in the basis 

f{Xs). 

Similarly, we can define the BI grid of the second type j/g by the two conditions 

-y2s=y2s+l, -1 - 2/2s = y2s-l, S = 0,±1,±2, . . .. (7.9) 

with the general solution 

{b — $; if s is even, 
—0 + II s IS odd, 

where b is an arbitrary real parameter. 
Equivalently, 

= -1/4+ (-1)^(1/4 + 6-5/2). (7.11) 
We then have that the operator H given by (7.4) is 3-diagonal as in (7.5) with 

^ ^ ( MV'^) if s is even, 

\B{y,) if sis odd, ^ ' 
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r Biy,) if s is even, 

if sis odd, ^ 

Vs = C{ys). (7.14) 
In particular, the BI difference equation (2.11) for the BI polynomials can be presented in the form 

£,sPn{Xs+i) + T]sPn{Xs) + CsPn{Xs-i) = XnPn{Xa), (7.15) 

where the coefficients ^sjVsjCs are given by (7.6)-(7.8) (or by (7.12)-(7.14)) with 

A{x) = -F{x), B{x) = G{x), C\x) = F{x) - G{x). 

In turn, property (7.15) means that the BI polynomials possess the Leonard duality property [3], [29]: they 
satisfy both a 3-term recurrence relation with respect to n and a 3-term difference equation with respect to the 
argument index s. Bannai and Ito found all orthogonal polynomials satisfying the Leonard duality property 
[3]. The corresponding classification, referred to as the Bannai- Ito theorem, is a generalization of the Leonard 
theorem [21] that deals only with polynomials orthogonal on a finite set of points. We have thus derived the 
Leonard duality property of the BI polynomials starting from the Dunkl shift difference equation (2.11). 
Note that the Bannai-Ito grid satisfies the linear equation 

Xs+i+Xs-i+2xs + l = (7.16) 

which is a special case of the generic equation 

X s+i + Xs-i - {q + q~^)xs = cxmst (7-17) 

characterizing the Askey- Wilson grids [14], [24]. As expected, the Bannai-Ito grid corresponds to the case 
q = —1 oi the Askey- Wilson grid. 

Consider now the special case when the operator H is symmetric H* = H, with real functions A{x), B{x), C{x). 
This is equivalent to the conditions 

A{-x)=A{x), B{-x -I) = Bix). (7.18) 

It is easily verified that in this case the 3-diagonal matrix in (7.5) becomes symmetric as well, i.e. = Cs+ij or 
in details: 

Hf{Xs) = Cs+lf{Xs+l) + VsfiXs) + Csf{Xs-l). (7.19) 

This result does not depend on the type of the BI grid. We already know that the operator M = ip{x)L 
is symmetric M* = M. Assume further that the function ip{x) is positive, i.e (fi{x) > for some interval 
xo < X < xi. We can then introduce the operator 

L = (p^/^{x)L^-^/^{x). (7.20) 

By condition (6.1), this operator will be symmetric L* = L on the interval [a;o,a;i]. 
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Assume additionally that the corresponding 3-diagonal matrix is finite, say it has dimension {N+1) x {N+1). 
This is equivalent to the conditions 

Co = Civ+i = (7.21) 
and Cs > 0) s = 1,2, . . . ,N. Let qn{s) be a set of eigenvectors of the s-diagonal matrix L: 

Cs+iqn{s + i) + VsQn{s)+Csqnis) = Kqnis), s,n = 0,1,2, N. (7.22) 

The operator L is a Hermitian matrix in the basis qn{s). By elementary linear algebra, it is well known that 
the coordinates qn (s) (under appropriate normalization) satisfy the orthonormality property 

N 

^gn(s)9m(s) = Snm- (7.23) 
s=0 

Prom the definition (7.20) of the operator L, it then follows that 

qn{s)=K-''fi'/^{Xs)Pn{Xs) (7.24) 

with some normalization factor k„. 

We thus arrive on the one hand, at the orthogonality relation for the BI polynomials: 

JV 

^ WsPn{Xs)Pm{Xs) = 5nm, (7.25) 

where the discrete weights are defined by the formula 

w, = ip{xs). (7.26) 

By construction, the weights are positive > as should be for positive definite orthogonal polynomials Pn{x). 

On the other hand, we have formula (4.44) giving the orthogonality relation of the BI polynomials in the 
finite case, where Xs are simple roots of the polynomial Pn+i{x) = {x — xo){x — xi) ... (a; — xn)- This means 
that the roots Xg can be parametrized in terms of the BI grid. 

Explicitly, the roots Xs of the polynomial Pn+i{x) can be found from formula (5.21) which expresses the BI 
polynomials in terms of the Wilson polynomials. Omitting technical details, we present the results. 

When A'' is even and condition (4.46) is assumed then 

a;, = -l/4+(-l)^(p2 + l/4 + s/2), s = 0, 1, 2, . . . , TV, (7.27) 

i.e. the roots for the case of even N belong to the BI grid of the first type (7.3) with a = p2- 
When N is odd and condition n + r2 = {N + l)/2 is assumed then 

X, = -l/4+(-l)«(-l/4 + ri -s/2), s = 0,1,2,..., N, (7.28) 
i.e. the roots for the case of odd N belong to the BI grid of the second type (7.11) with 6 = n — 1/2. 
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So for even AT = 2, 4, 6, . . . the orthogonality relation reads as in (7.25) with the spectral points Xg given by 
(7.27) and the discrete weights expressed as in (7.26), and where ip{x) can be taken in the form (6.13). 

For odd N = 1, 3, 5, . . . the spectral points Xs are given by (7.27) and the discrete weights expressed as in 
(7.26), where ip{x) can now be taken of the form (6.14) (this form is preferable for the case of odd N in order 
to avoid singularities when x = Xg)- 

It is easily verified that in all cases the function (7o{xs) and (T2{xs) are constant sign-alternating sequence, 

e.g. 

cro{xs) = {-lymnst, s = 0,1,2, ... ,N 

Note finally that the set of spectral points Xs is the union of two discrete subsets X2s and X2s+i corresponding 
to even and odd values of s. These two subsets never overlap if the positivity conditions form the BI polynomials 
are fulfilled. 



8. Bannai-Ito and complementary Bannai-Ito polynomials as limits 
g — > — 1 of the Askey- Wilson polynomials 

The Bannai-Ito polynomials were introduced in [3] through a direct limit process from the g-Racah polynomials 
as g —1. 

We show in this section how the Bannai-Ito polynomials can be obtained from the Askey- Wilson (AW) 
polynomials as q — )• —1. We also indicate how the complementary Bannai-Ito polynomials can equally be 
derived from the AW-polynomials in the same limit. 

Consider the Askey- Wilson polynomials ii„(x(z); a, h, c, d) [2], [17] 

where 

'^1^'' =^ a^^'iab, ac, ad; q)n- 

The polynomials Rnix{z): a,b, c, d) depend on the argument x = {z + z~^)/2 and on 4 complex parameters 
a,b,c,d. They satisfy the recurrence relation [17] 

Ar,Rn+i{x) + {a + a-^ - A„ - C„)i?„(a;) + C„i?„_i(a;) = {z + z-'^)Rn{x) = 2xRn{x) (8.2) 

with coefficients 

(1 - a5g")(l - acg")(l - a(ig")(l - abcdq''-^) 



a{l - a6cdq'2"-i)(l - abcdq^'<^) 



^ a(l - g")(l - 6cg"-i)(l - Mg"-i)(l - cdg"-^) 

(1 - a6c(ig2"-i)(l - o6cd(72"-2) " ^"-^^ 

In order to recover the Bl-polynomials we choose the following parametrization 

q = -e^ a = -ie'", b = -ie^^ , c = ie^^ , d = ie^\ (8.4) 
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where a, 13, 7, S are real parameters. In what follows these parameters are more conveniently expressed in terms 
of the 4 parameters ri,2,pi,2: 

a = 2pi + l/2, /3 = 2p2 + l/2, 7 = -2r2 + 1/2, 5 = -2ri + 1/2. (8.5) 

Moreover, we take 

z = ie-'^'y (8.6) 

with y a real parameter. 

It can be checked that the limit g — >■ — 1 of the Askey- Wilson polynomials 

lim Rn{x{z)) = R^-^\y) (8.7) 

g-)— 1 

does exist, and that Rii^^\y) is a polynomial of degree n in the argument y. 

Dividing the recurrence relation (8.2) by 1+q and taking afterwards the limit e — > (which is equivalent to 
the limit g — )• —1), we obtain a recurrence relation for the limiting polynomials R^^\y): 

4"'^^i;V(2/) + (1/4 + P1- 4"'^ - c(-i))i?(-i)(2/) + cf^Rt'hy) = yRi-'Hy), (8.8) 

where the coefficients A^^\ Ci~^^ have expressions coinciding with (4.41) and (4.42). The polynomials i?i~^^(i/) 
are thus identified with the Bannai-Ito polynomials (up to a shift of the argument): 

Ri-'\y)=P„{y-l/A;ri,r2,pi,P2). (8.9) 

R^^\y) in (8.9) denotes the monic version of the polynomials i?i~^'(t/), while P„(y; ri, r2, pi, P2) stands for the 
monic Bannai-Ito polynomials defined in (4.40). The parameters ri,r2, Pi, P2 for the Bannai-Ito polynomials 
stem under the limiting procedure (8.5), from the corresponding 4 parameters of the Askey- Wilson polynomials. 

For the complementary Bannai-polynomials, we use the same procedure but with a slightly different parametriza- 
tion: 

q = -e^ a = ze^^^^^+^Z^), b = -ie<P-+^/^\ c = ie'^-^^-+^/^\ d = ie'^-^^'+^/^\ z = ie-^^y. (8.10) 

The main difference between (8.4) and (8.10) is a change of sign in a and the shift pi — >■ pi + 1/2. 

It is directly verified that in the limit £ — > 0, we get the recurrence relation (5.12) that defines the CBI 
polynomials and therefore that 

\im Rn{z; a, b,c,d) = Wn{y - 1/4; n, r2, pi, P2). (8.11) 

Thus the BI and CBI polynomials are obtained through very similar q ^ —1 limits of the Askey- Wilson 
polynomials. 
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9. The Bannai-Ito Dunkl shift operator as a hmiting form of the 
Askey- Wilson difference operator 

In this section, we describe how the Bannai-Ito Dunkl shift operator appears in the limiting process q ^ —1 
from the difference Askey- Wilson operator. 

The Askey- Wilson polynomials Rn{x{z)) satisfy the difference equation [17] 

n{z)Rn{x{zq)) + n{z-^)Rn{x{zq-^)) - {n{z) + Vl{z-^))Rn{x{z)) = AnRn{x{z)), U = 0, 1, 2, . . . , (9.1) 

where 

A„ = (g-" - 1)(1 - abcdq"--'^) 

is the eigenvalue and 



"(^) = T. • (9-2) 



(1 - az){l - bz){l - cz)(l - dz) 

(l-z2)(l-gz2) 

The linear operator in the Ihs of (9.1) (acting on the variable z) is the Askey- Wilson difference operator [2], 
[17]. 

Consider the limiting form of the difference equation (9.1) when q — ^ —1. As in the previous section, we 
choose the parametrization (8.5) and (8.6). Wc already showed that the Askcy-Wilson polynomials Rn{x{z)) 
become the Bannai-Ito polynomials i?l ^\y)- The operation Rn{z) Rn{zq) is reduced to the operation 
Rri^^\y) Rri^^\—y + 1/2) while the operation i?„(z) — >■ Rn{zq^^) becomes the operation Rii^^\y) — >■ 
Ri-'\-y- 1/2). 

It is then straightforward to verify that in the limit £ ^ we get the equation 

^^{y)Ri-'\-y + 1/2) + My)Ri-'\-y - 1/2) - ($i(y) + ^2{y))Ri-'\y) = KR^n^Ky), (9.3) 

where 

e^0 4(l-Fg) Ay -I 

My) = lim = ^iy-n + miy-r. + m 

e^o4(l + g) 4y+l ^ ' 

and that the eigenvalue 

Xn = lim , (9.5) 

s^oA{l + q) ^ ' 

has the expression (2.13). 

This gives the Bannai-Ito polynomials R^^\y) as eigenfunctions of some Dunkl shift operator. In operator 
form this eigenvalue equation can be presented as follows: 

HRi-^\y) = X^Ri-^\y). (9.6) 

H in (9.6), acts on the space of functions f{y) of argument y according to 

H = My)T-^'^R + My)T^/^R - ($i(y) + ^y)), (9.7) 
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and stands for the shift operator 

T^f{y) = f{y + h). 

The Dunkl shift operator H contains (apart from the identity operator) the operators T^/^jT^^/^ and R. 

We can now reduce the operator (9.7) to the "canonical" form which only involves the operators T+ and R. 
It goes as follows. The similarity transformation 

H = T^HT-^ 

with an arbitrary real /i is a unitary transformation of the operator H (recall that by T"*" and T~ we mean the 
operators and respectively). Under this transformation H goes into 

fl- = $i(y + h)T^^-''/^R + ^2{y + h)T^^+^^^R - ^{y + h)-^i{y + h). (9.8) 

For h = 1/4, we have 

H = ^i{y + l/4)i? + $2(y + l/4)T+i? - $i(y + 1/4) - $2(y + 1/4), (9.9) 

where 

$,(,+ l/4) = -fc^ifc^ = -F(,), ^,iy + l/4) = Jy^^l±l^i^f:^±ll^ = Giy). (9.10) 

The functions G{x),F{x) in (9.10) are precisely those of (2.12). We have thus reduced the operator H to our 
initial "canonical" Dunkl shift operator L defined in (2.4). 

Note that there is one more "canonical" form of the Dunkl shift operator where only T~ and R appear. 
This other form is obtained by choosing ft = —1/4, so that 

H = ^i{y- l/A)T-R + $2(t/ - l/4)i? - $i(y - 1/4) - $2(2/ - 1/4). (9.11) 

We see that the Dunkl shift operator can be presented in 3 "canonical" form: the "+ form" given by (9.8), the 
"— form" given by (9.11) and the "symmetric" form given by (9.7). Clearly, all these forms are equivalent and 
correspond to shifts y ^ y + h m. the argument of the polynomials. 

As was already mentioned, the complementary BI polynomials do not satisfy a difference equation of Dunkl 
shift type. Indeed, under the choice (8.10), the corresponding Askey- Wilson difference operator does not survive 
in the limit g — >■ — 1. 

10. Symmetric BI polynomials and continuous dual Hahn polyno- 
mials 

We already saw that for generic BI polynomials with real parameters ri,r2, pi, P2 it is impossible to obtain 
positive definite polynomials for all n = 0, 1, 2, 
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There is however, a special case for which it is possible to construct such positive definite polynomials after 
the following change of argument: x ix. This case arises for symmetric BI polynomials. 
Recall that polynomials Sn{x) are said to be the symmetric if they satisfy the condition 

= (-1)" 5„(x). (10.1) 

Orthogonal polynomials Pn{x) obeying a recurrence relation of the form 

Pn+l{x) + (x) + U„P„_i(x) = xPn{x) 

are symmetric iff 6„ = [6] (i.e. iff the diagonal recurrence coefficient 6„ vanishes for all n = 0, 1, 2, . . .). 

If 6„ = 6 = const, n = 0, 1, 2, ... it is then the "shifted" polynomials Pn{x) = Pni'x + h) that are symmetric. 

Let us determine when the (shifted) BI polynomials B!i^^\x) (defined by (8.9)) are symmetric. This is 
equivalent to determining when the condition 6„ = — 1/4, n = 0, 1, 2, ... for the diagonal recurrence coefficient 
of the BI polynomials, is satisfied. It is easy to see that this condition holds iff 

either 

(i) pi = -n, p2 = -r2; 
or 

(ii) pi = -r2, P2 = -n. 
In both cases we have 



and 



6„ = -l/4, n = 1,2,3,... (10.2) 

^ J_ n(n-4(ri+r.)) ^ if 71 evCn 
"»-|_ (n-4..)(n-4..) ^ if n odd' ^^^'^^ 

The corresponding symmetric polynomials Sn{x) = Ri^^\x) = Pn{x — 1/4) are expressed in terms of dual 
Hahn polynomials. In order to obtain positive definite polynomials, we need to perform the change of variable 
X — >■ ix. This leads to symmetric polynomials with positive coefficients «„. Indeed, for the new polynomials 
Sn{x) = i~^Sn{ix), we have the 3-term recurrence relation 

Sn+l (^) + ) = xSn{x) (10.4) 

with coefficients Un = i^Un — —Un and thus 



— " . if n even 

if n odd 



(n-4,.).,-4,.) (10-5) 



I 16 

It is seen that under the restrictions 

n < 1/4, ra < 1/4, 

one has u„ > for all n = 1,2,3, This means that the polynomials Sn{x) are orthogonal on the real axis 

with a positive measure. 
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Comparing the recurrence coefBcients (10.5) with those for the continuous dual Hahn polynomials [17] and 
using Chihara's standard procedure to relate symmetric and non-symmetric polynomials [6], we can derive an 
explicit expression for the polynomials Sn{x): 

([)•) I —n,a + 2ix,a — 2ix A , , (t\ ^ f —n,a + 2ix,a — 2ix 

CA"' o /<o '1 ,Sn_ 1 1 I Tl =: fcA ' ^ I' ' 



(10.6) 

where 

a = -2ri + 1/2, b= -2r2 + 1/2 

and the normalization coefficients are 

4°) = (-l)"2-2" («)„ (a + b)n, = (-l)"2-2" (a + 1)„ (a + 6)„. 

Recall that the (unnormalized) continuous dual Hahn polynomials Cn(a;^; a, b, c) depend on 3 parameters a, b, c 
and are expressed as [17] 

^ / 2 , \ ^ f -n,a + ix,a - ix \ /,„„x 
Cn{x^;a,b,c) = 3F2[ ' , ' ;1 • 10.7 
\ a + o,a + c J 

We thus see that 

S2n{x) - 4"^ Cj4x^;a,b,0), S^n+iix) - a;C„(4a;2; a, &, 1). 

From the expression of the weight function of the continuous dual Hahn polynomials [17], we reconstruct the 
weight function of the symmetric BI polynomials 

, , T{a + 2ix)T{b + 2ix) ^ 
^("^ = r(l/2 + 2ix) 

for which they are orthogonal on the whole real line 
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f 

J — c 



'w{x)Sn{x)Sm{x)dx = 0, Tl ^ m. 



The difference equation for these polynomials follows immediately from the corresponding difference equation 
(9.3) for the BI polynomials after replacing y ^ iy: 

My)Sn{-y - + My)Sni-y + 1/2) - ($1(2/) + <^>2{y))Sn{y) = A„ Sn{y), (10.9) 

where 

,^,^,_^(«^, .,„.. (»-^/2K.- .v.) 

and the eigenvalue is 

if n is even 

2(,,+,,)_ln+li if n is odd. ^'^-'^^ 
Consider the special case when one of the parameters ri or r2 vanishes, say r2 = 0. The recurrence coefficients 
Un then become 

.^^n(n_4ri)^ n= 1,2,3,... (10.11) 
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These correspond to the recurrence coefficients of the symmetric Meixner-Pollaczek polynomials [17] 

The corresponding Dunkl shift eigenvalue equation is 

~{2y + ia)Sn{-y - + ^-{2y - ia)Sn{-y + i/2) + | 5„(y) = A„ Sn{y). (10.13) 

Taking into account the property Sn{—y) = {—l)^Sn{y) of the symmetric polynomials, equation (10.13) can be 
cast in the form 

--{2y + ia)Sn{y + 1/2) + -{2y - ia)Sn{y - %I2) = ^-5„(y) (10.14) 

which coincides with the ordinary difference equation of the Meixner-Pollaczek polynomials on a uniform imag- 
inary grid [17]. 

The weight function in this case is obtained from the weight function (10.8) by putting h = 1/2: 

w{x) = \^[a^2ix)f , 

which coincides with the well known expression for the weight function of the symmetric Meixner-Pollaczek 
polynomials [17]. 

The symmetric Meixner-Pollaczek polynomials therefore belong, as a special case, to the class of BI poly- 
nomials. The eigenvalue equation for these polynomials can be presented either in the form (10.13) as an 
eigenvalue equation for a Dunkl shift difference operator of first order, or in the form (10.14) which is a stan- 
dard Sturm-Liouville difference equation of the second order. 



11. The big and little -1 Jacobi polynomials as limit cases 

Changing the variable to x = y/h with h a real parameter, we can rewrite the BI operator in the form 

Lfiy) = fc^l||^(/(,) - fi-y)) + ^'-'''^ll^y%'''^'^'\ fi-y - h) - m). (11.1) 
Choose now the following parametrization 

ri = ai/h, r2 = a2/h, pi = ai/h + bi, p2 = a2/h + 62, (11-2) 

where the real parameters ai,bi,i = 1,2 do not depend on h. Taking the limit /i — >■ 0, we then arrive at the 
operator 

where 

F{y) = (26i + 262 + l)y^ - 2(0163 + a2bi)y - aia2. 
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The operator (11.3) is a first order differential operator of Dunkl type. It was shown in [34] that the operator 
Lo as given in (11.3) with ai,a2,bi and 62 real, represents the most general first order differential operator of 
Dunkl type that has orthogonal polynomials as eigenfunctions. In general the operator Lq has 4 free parameters. 
However, if 02 7^ 0, by scaling the argument y — > 7?/ one can set 02 = 1 and only 3 independent real parameters 
01,61,62 remain. 

Using a slightly different parametrization, we thus have [33] 

Lo = go{y){R-l)+gi{y)dyR, (11.4) 

where 

ia + P + l)y^Hca-P)y + c ^ ^ 2(,-l)(, + c) 

y y 

The 3 independent real parameters of the operator Lq are now a, j3 and c. The polynomial eigensolutions 
Pjiiy), n = 0, 1, 2, . . . of io are expressed in terms of the big -1 Jacobi polynomials (see [33] for details). 
If 02 =0, the operator Lq can be reduced to 

Lo = 2(1 - y)dyR + (a + /3 + 1 - ay-^){l - R) (11.6) 

with only two independent real parameters a, 13. The polynomial eigensolutions Pn{y),n = 0,1,2,... of the 
operator (11.6) coincide with the little -1 Jacobi polynomials [32]. 

In [33] it was shown how the recurrence coefficients of the big -1 Jacobi polynomials could be recovered 
through a special large N limit, from those of the Bannai-Ito polynomials. Clearly, the eigenvalue problem 
for the big -1 Jacobi polynomials can also be obtained from the corresponding eigenvalue problem for the BI 
polynomials. For g — > 1 this corresponds to the transition from the q-Racah to the big q-Jacobi polynomials 
[20]. 



12. Conclusion 

The study of the polynomial eigenfunctions of first order Dunkl shift operator has thus led synthetically to a 
thorough characterization of the Bannai-Ito polynomials: weight function, structure and recurrence coefficients 
etc. It has brought to light a relation between certain Jordan algebras, their representations and orthogonal 
polynomials. It made natural the introduction of the complementary Bannai-Ito polynomials, their expression 
in terms of Wilson polynomials and their connection to the BI polynomials. It further allowed to see that the 
symmetric BI polynomials are related to the dual Hahn polynomials and to study finally, various limiting eases. 

Note that polynomial eigenfunctions of certain linear operators with refiections were previously considered in 
the literature (see for instance [19]); these solutions however exhibit non-standard orthogonality properties. In 
contrast, the polynomial eigensolutions that have arisen in our treatment, are clearly orthogonal in the ordinary 
sense. 
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